Closed-shell properties of ^''O with ab initio coupled-cluster theory 
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We present an ah initio calculation of spectroscopic factors for neutron and proton removal from 
^■'O using the coupled-cluster method and a state-of-the-art chiral nucleon-nucleon interaction at 
next-to-next-to-next-to-leading order. In order to account for the coupling to the scattering contin- 
uum we use a Berggren single-particle basis that treats bound, resonant, and continuum states on 
an equal footing. We report neutron removal spectroscopic factors for the '^^O states = 1/2^, 
5/2"*", 3/2~ and l/2~, and proton removal spectroscopic factors for the '^'^N states l/2~ and 3/2~. 
Our calculations support the accumulated experimental evidence that ^''O is a closed-shell nucleus. 

PACS numbers: 21.10.Jx, 21.60.De, 21.10.Pc, 31.15.bw, 24.10.Cn 



Introduction The study of nuclei far from stability 
is a leading direction in nuclear physics, experimentally 
and theoretically. It represents a considerable intellectual 
challenge to our understanding of the stability of matter 
itself, with potential implications for the synthesis of ele- 
ments. An important aspect of this research direction is 
to understand how magic numbers and shells appear and 
evolve with increasing numbers of neutrons or protons. 
The structure and properties of barely stable nuclei at 
the limits of stability, have been demonstrated to deviate 
dramatically from the established patterns for ordinary 
and stable nuclei, see Ref. [l| and references therein for a 
recent review. One of the striking features of nuclei close 
to the drip line is the adjustement of shell gaps, giving 
rise to different magic numbers 0] . The way shell closures 
and single-particle properties evolve as functions of the 
number of nucleons forms one of the greatest challenges 
to our understanding of the basic features of nuclei, and 
thereby the stability of matter. 

The chain of oxygen isotopes up to ^^0 is particularly 
interesting since these are the heaviest nuclei for which 
the drip line is well established. Two out of four stable 
even-even isotopes exhibit a doubly magic nature, namely 
(Z=8,7V=14) and ^^O (Z=8, N=l&). Several recent 
experiments bring evidence that ^"'O is the last sta- 
ble oxygen isotope. This is remarkable, in particular if 
one considers the fact that the addition of a single pro- 
ton on top of the Z — % closed shell brings the drip line 
of the fluorine isotopes to '^^F. Recent measurements [3] 
also suggest that ^^O has a spherical neutron configura- 
tion. The isotopes 25-28q g^j.g g^jj believed to be unstable 
towards neutron emission, even though ^^O is a doubly 
magic nucleus within the standard shell-model picture. 
This indicates that the magic number at the neutron drip 
line for the oxygen isotopes is not at iV = 20 but rather 
at TV = 16. 

Although S2ectroscopic factors are not observable 
quantities [1, [13] , they can be used to address shell clo- 
sure properties within the context of a given model. Ex- 
perimentally, spectroscopic factors are defined as the ra- 
tio of the observed reaction rate with respect to the same 



rate, calculated with a particular reaction model and as- 
suming a full occupation of the relevant single-particle 
states. The spectroscopic factors are then interpreted as 
the occupancy of specific single-particle states. Theoret- 
ically however, spectroscopic factors measure what frac- 
tion of the full wave function corresponds to the product 
of a correlated state (often chosen to be a given closed- 
shell core) and an independent single-particle or single- 
hole state. Large deviations from the values predicted by 
an independent-particle model point to a strongly corre- 
lated system. 

Kanungo et al reported measurements of one- 
neutron removal from ^^Q, and an extraction of spec- 
troscopic factors. They used an eikonal reaction model 
and Woods-Saxon overlap functions, that were calculated 
with a well depth adjusted to reproduce the one-neutron 
separation energies. The results were compared with 
shell model predictions using the fitted sc?-shell USDB 
interaction jj] and the SDPF-M interaction Q . The theo- 
retical calculations corroborate the large s-wave probabil- 
ity found in the experimental analysis, implying thereby 
that is indeed a doubly magic nucleus. 

The aim of this work is to add further theoretical ev- 
idence and support to these claims. We present spec- 
troscopic factors for neutron and proton removal from 
2*0 as predicted by the ah initio coupled-cluster method 
with the chiral nucleon-nucleon interaction by Entem 
and Machleidt JTj at next-to-next-to-next-to-leading or- 
der (N'^LO). Our calculations are performed in a large 
single-particle basis which includes bound and continuum 
single-particle states. The above mentioned theoretical 
calculations of Refs. 0, Q involve only fitted effective in- 
teractions tailored to small shell-model spaces, such as 
the sd or the sd — pf shells only. 

The virtue of ab initio methods applied to nuclear 
physics is to reduce the model dependence of computed 
results. By means of a recipe for systematic improve- 
ments, one can distinguish between parameters of tech- 
nical and physical character. Whereas the result may 
be expected to depend on physical parameters, it should 
be insensitive to the technical parameters as systematic 
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improvements are included. Ultimately, a converged ab- 
initio result may provide a rigorous test of the nuclear 
interaction model and the corresponding physical param- 
eters. Since our single-particle basis contains continuum 
states, and therefore represents the correct asymptotical 
behaviour, we are able to generate radial overlap func- 
tions for drip-line nuclei. This work paves therefore the 
way for neutron-knockout reactions with fully ab initio 
structure information. 

After these introductory remarks, we briefly expose our 
calculational formalism as well as the choice of Hamilto- 
nian and basis in the next sections. Following that, our 
results are presented. Conclusions and perspectives are 
drawn in the final section. 

Method The aim of this section is to give a short 
overview of the step s in our calculations. The coupled- 
cluster method I12h22| | and the application to spectro- 
scopic factors |23[ have been presented in great detail 
elsewhere. Here we give only a brief summary of the 
concepts that enter the calculations leading to the results 
presented in this article. 

The spectroscopic factor (SF) is the norm of the over- 
lap function, 

SU(.lj) = \0Uilj;r)\\ (1) 
Oi_,ilr, r) = ${A-l\\ hra, II A) Ki,{r) . (2) 

n 

Here, 0^_i{lj\r) is the radial overlap function of the 
many-body wavefunctions for the two independent sys- 
tems with A and A — 1 particles respectively. The double 
bar denotes a reduced matrix element, and the integral- 
sum over n represents both the sum over the discrete 
spectrum and an integral over the corresponding contin- 
uum part of the spectrum. The annihilation operator 
<inijm transforms like a spherical tensor of rank j and 
projection to. The radial single particle basis function is 
given by the term (pnij (r) , where I and j denote the sin- 
gle particle orbital and angular momentum, respectively, 
and n is the nodal quantum number. The isospin quan- 
tum number has been suppressed. We emphasize that 
the overlap function, and hence also its norm, is defined 
microscopically and independently of the single particle 
basis. It is uniquely determined by the many-body wave 
functions \ A) and |^ — 1). The quality of the SF estimate 
is thus limited by the quality of the pertinent many-body 
wave functions. 

Our calculation of spectroscopic factors follows the 
recipe detailed in Ref. [1^. The important difference 
between this work and Ref. [1^, is that all terms con- 
tributing to the spectroscopic factors have now been ex- 
pressed in terms of reduced matrix elements in an angu- 
lar momentum coupled basis. This allows us to handle 
a much larger set of single-particle states as discussed in 
Refs. 0,1211 . 

We use the coupled-cluster ansatz, — exp (T) |0o} 
for the ground state of the closed-shell nucleus ^''O. The 
reference state, |0o}: is an antisymmetric product state 
for all A nucleons. The cluster operator T introduces 



correlations as a linear combination of particle-hole exci- 
tations T — Ti -\- T2 -\- ■ ■ ■ -\- Ta^ where T„ represents an 
ri-particle-n-hole excitation operator. For the coupled- 
cluster singles and doubles approximation (CCSD) em- 
ployed in this work, T is truncated at the level of double 
excitations, T = Ti + T2. The coupled-cluster solution 
for ^■^O is obtained as a set of amplitudes that defines T. 

Due to the non-hermiticity of the standard coupled- 
cluster ansatz, we need both the left eigenvectors and 
the right eigenvectors. These are determined via 
the cquation-of-motion coupled-cluster (EOM-CC) ap- 
proach as \A) « \R^{Ja)) = exp{T)R^{JA)\(t)o) and 
{A\ « {L^{Ja)\ = {<j>o\L^{JA)exp{^T). The opera- 
tors R^{Ja) and L^{Ja) produce linear combinations 
of particle-hole excited states when acting to the right 
and left, respectively. In the spherical form of the EOM- 
CC approach, the operators have well defined angular 
momentum by construction, as indicated by Ja, which 
stands for the angular momentum considered. See for 
example Ref. [2l| . If the A-body system is in its ground 
state, the right EOM-CC wave function is identical to 
the coupled-cluster ground state. 

Solutions for the A — 1-body system, ^'^O and ^^^N, are 
obtained with particle-removal equations-of-motion (PR- 
EOM-CCSD), where we use the CCSD ground state so- 
lution of ^^O as the reference state in order to determine 
the corresponding left and right eigenvectors as |^ — 1) ~ 
\R^-HJa-i)) = exp iT)R^'HJA-i) \M and {A - 1\ ^ 
{L^-\Ja-i)\ = {(l)o\L^-\JA-i)exp{-T). In actual 
calculations, the EOM-CC wave functions are obtained 
by determining the operators R^{Ja) and L^{Ja) as 
eigenvectors of the similarity transformed Hamiltonian, 
H — exp (— T)iJ exp (T). The transformed Hamiltonian 
is non-hermitian, implying that the left eigenstates must 
be determined independently. We refer the reader to 
Refs. [13, 1131 for details about the equation-of-motion 
approach combined with coupled-cluster theory. 

Finally, we can approximate the spectroscopic factor, 
Eq. ([TJ in the spherical coupled-cluster formalism as 

Si-iilj) = t {L^-'iJA^i)\\K^\\R^iJA)) 

x{Rfi-\jA-i)\\h^A\Lt{JA))\ (3) 

where we have used the similarity transformed spherical 
annihilation operator as 

anijm = exp {-T)dnij„i exp (T). 

Closed expressions for the similarity transformed opera- 
tors are given in Ref. [23]. The labels fi and v are included 
to distinguish excited states of the (PR-)EOM-CC solu- 
tions. In the spherical formulation of EOM-CCSD, the 
solutions are spherical tensors 21], and the spectroscopic 
factor depends on the rank, but not on the projection of 
the EOM-CCSD states. In order to derive the coupled 
expressions, a Racah algebra module was developed for 
the open source computer algebra system SymPy [13]. 
More details about these calculations can be found in 
Ref. [13]. 
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Hamiltonian and Basis We use an intrinsic A-nucleon 
Hamiltonian H = T — + V where T is the ki- 
netic energy, Tcm is the kinetic energy of the center-of- 
mass coordinate, and V is the two-body nucleon-nucleon 
interaction. Coupled-cluster calculations starting from 
this Hamiltonian have been shown to generate solutions 
that are separable into a gaussian center-of-mass wave 
function and an intrinsic wave function, see for example 
Refs. dH, m for further details. 

The nucleon-nucleon interaction we use is the chiral 
N'^LO interaction model of Entem and Machleidt 
Although the interaction has a cut-off A — 500 MeV, it 
still contains high momentum modes, and one typically 
needs model spaces which comprise about 20 major oscil- 
lator shells in order to reach convergence for the ground 
states of selected oxygen and calcium isotopes, see for 
example Refs. [Ulll^ for a discussion. However, exploit- 
ing the spherical symmetry of the interaction and our 
coupled-cluster formalism, we can use model spaces that 
are large enough so that there is no need for a subsequent 
renormalization of the interaction. 

Short-range properties are not modelled explicitly, but 
are instead represented by contact terms in the interac- 
tion. The contact terms are subject to the fitting proce- 
dure, meaning that the short range dynamics will be sen- 
sitive to the momentum cut-off. Since we are neglecting 
many-body forces like three-body forces or more compli- 
cated terms, our results will in general depend, less or 
more, on the chosen cut-off of the nucleon-nucleon inter- 
action model. Quantities that are sensitive to the short- 
range part of the wavefunctions, such as spectroscopic 
factors, may depend strongly on the chosen cut-off. This 
is why even ah initio calculations of spectroscopic factors 
must be considered as model dependent. Our results may 
be fully converged in terms of a given Hamiltonian and its 
parameters at a given level of many-body physics, how- 
ever, employing another nucleon-nucleon interaction may 
lead to slightly different results since many-body terms 
beyond those represented by a two-body interaction can 
be very important. This is discussed in detail in Ref. [25j . 

We use a Gamow-Hartree-Fock (GHF) solution for the 
reference state, as detailed in for example Ref. 20] . These 
Hartree-Fock solutions were built from the standard har- 
monic oscillator (HO) basis combined with Woods-Saxon 
(WS) single particle states for selected partial waves in 
order to properly reproduce effects of the continuum. 
The role of the continuum is expected to be important 
close to the drip line, as seen in Refs. [13, [26l - [28|. For 
this purpose we use a Berggren representation '29'] for 
the neutron Si/2, (^3/2; a-nd partial waves. This rep- 
resentation generalizes the standard completeness rela- 
tion to the complex energy plane. In the Berggren basis, 
bound, resonant, and non- resonant continuum states are 
treated on an equal footing. The Berggren ensemble has 
been successfully used within the Gamow shell model, see 
for example Ref. [27| for a recent review, and in ah ini- 
tio coupled-cluster calculations of energies and lifetimes 
in Refs. [1^ [l^]. The Berggren basis is constructed by 



diagonalizing a one-body Hamiltonian with a spherical 
Woods-Saxon potential in a spherical plane-wave basis 
defined on a discretized contour in the complex momen- 
tum plane. We employ a total of 30 Gauss-Legendre mesh 
points along the contour for each of the Si/2, ^3/21 and 
d^l2 neutron partial waves. With 30 discretized single- 
particle states our results for the above single-particle 
states become independent of the choice of contour. For 
the choice of interaction that we use, the 1/2+ and 5/2+ 
states are fairly well bound with respect to ^^O [2ll[3l|. 
therefore it is sufficient to use a contour along the real en- 
ergy axis. For all other partial waves, the basis functions 
are those of the spherical harmonic oscillator. 

Results Figured] shows the spectroscopic factors for re- 
moving a neutron in the Si/2 and d^/2 partial waves of 
^^O as function of the harmonic oscillator frequency Sw. 
The ground state 1/2+ and excited 5/2+ state in ^^O 
were calculated within the PR-EOM-CCSD approxima- 
tion starting from the GHF basis with 30 mesh points for 
each of the si/2, c?3/2 and ^5/2 neutron partial waves and 
17 major oscillator shells for the protons and remaining 
neutron partial waves. The spectroscopic factors are well 
converged with respect to the model space size. To inves- 
tigate the role of continuum on the spectroscopic factors, 
we compare with a calculation done in a hartree-fock ba- 
sis built from a harmonic oscillator basis of 17 major 
shells (OHF). 

We find that the effect of the continuum is small. This 
is expected, since our calculations of the 1/2+ and 5/2+ 
single particle states in ^"'O result in well bound states 
with respect to_the neutron emission threshold, see for 
exmaple Ref. 
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for more details. We do however see 
a small reduction of the spectroscopic factors when the 
continuum is included. The reduction of spectroscopic 
factors will be enhanced for states close to a reaction 
channel threshold as discussed in Ref. tSSi] . Although the 
effect of the continuum on the spectroscopic factors is 
marginal in the cases we consider, the effect is crucial 
in order to obtain the correct asymptotic behaviour of 
the overlap functions for one-neutron removal. It is the 
asymptotic normalization coefficient, which is calculated 
from the tail of the radial overlap function that enters 
the exact reaction amplitudes. It is arguably the rele- 
vant quantity to calculate for reactions |34| . The results 
in Fig. [T] show also that the spectroscopic factors for the 
single neutron states close to the Fermi surface are close 
to one, indicating that a configuration consisting of a 
single-hole removal from ^'^O captures much of the struc- 
ture of the wave function for these states. The results for 
the neutron si/2 and d^/2 hole states lend thus support to 
the accumulated evidence that ^^O can be interpreted as 
a closed-shell nucleus. For the proton states close to the 
Fermi surface, this finding is only partly corroborated by 
the results shown in Fig. [21 In that figure we show our 
results for the spectroscopic factors for removing either a 
neutron or a proton in the P3/2 and pi/2 partial waves of 
^''O. From Fig.[5]we see that the 3/2~ and l/2~ states in 
neither ^■^O nor ^'^N can be clearly interpreted as simple 
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FIG. 1. (Color online) Normalized spectroscopic factors plot- 
ted against fujj for J'" = 1/2+ and 5/2+ one- neutron removal 
from ^■'O. The continuum states included in the Berggren- 
basis calculation (GHF), leads only to a small reduction com- 
pared with the harmonic oscillator values (OHF). The depen- 
dence on the oscillator spacing f'lw is very weak in both calcu- 
lations. The inclusion of the continuum structure also gives a 
small, but visible, improvement in the hio dependence. 
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for the 1/2+ state, and the theoretical results obtained 
with fitted shell model interactions. The USDB inter- 
action gives 1.810 for the neutron 1/2+ state and 5.665 
for the neutron 5/2+ state. Corresponding numbers for 
SDFP-M are 1.769 and 5.593, respectively [a]. 

As can be seen from Table HI our spacing between the 
two one-hole states is only 0.35MeV. This contradicts 
the shell model calculations with fitted interactions re- 
ported in Ref. Q. The energy spacing is 2.586MeV and 
2.593MeV for the SDFP-M interaction and the USDB 
interaction, respectively. Assuming a large spectroscopic 
factor for the 5/2+ state, a bound 5/2+-wave state at 
an energy consistent with the shell model predictions, 
should have been seen in the experiment of Kanungo et 
al f5] . The lack of such an observation was interpreted as 
a confirmation that this state is unbound. We speculate 
that missing three-nucleon forces could play an impor- 
tant role regarding the shell gap between the 1/2+ and 
5/2+ single particle states in ^^O. However, recent calcu- 
lations by Otsuka et al [3l| where three-body interactions 
were included as density-dependent corrected two-body 
interactions in a shell-model calculation, give results for 
^"^O which are very close to our spacing of 0.35 MeV. 
Similarly, a many-body perturbation theory calculation 
using the same Hamiltonian as here, results in a spac- 
ing of 0.39 MeV. Further investigations are thus needed 
in order to clarify the discrepancies between the results 
reported in for example Ref. Q, the present results and 
those of Otsuka et al [sij. Compared with the experi- 
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FIG. 2. (Color online) Normalized spectroscopic factors plot- 
ted against fioj for the negative parity states = l/2~ and 
3/2" in one-proton (vr) and one-neutron iy) removal from 
24o. 



one-hole states in O. However, the stronger dependence 
on the harmonic oscillator frequency huj indicates that we 
are missing many-body correlations beyond the 2h — \p 
level in our PR-EOM-CCSD computations of the 1/2" 
and 3/2- states. 

Our results for spectroscopic factors and energies for 
the J"" = l/2+,5/2+,3/2-, 1/2- states in and the 
J'^ = l/2-,3/2- states in ^^N are summarized in Ta- 
ble m We present the values that were calculated with a 
harmonic oscillator frequency fiw = 30 MeV, and the en- 
ergies are given relative to the ground state of ^"^O. The 
spectroscopic factors agree with both the extracted value 



TABLE I. Energies for states in '^O and N and correspond- 
ing spectroscopic factors (SF) for the removal of a particle 
from 2^0. The reported coupled-cluster results are at the 
level of singles and doubles (CCSD), and are calculated with 
fiw = 30MeV. Experimental values for ^^q are taken from 
Ref. [I and data for is taken from Ref. 

mental results in Table HI we note that the experimen- 
tal error bars are typically orders of magnitude larger 
than the dependence on technical parameters displayed 
in Figs. [T] and [2l Finally, we note that the 1/2- state in 
•^■^N is in very good agreement with the theoretical mass 
evaluations of Ref. ^] . 

However, there is still a considerable model dependence 
inherent in the calculation of spectroscopic factors. First 
of all, our calculations have been performed at the level 
of the singles and doubles approach, meaning that all cor- 
relations up to the level of two-particle-two-hole excita- 
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tions are included to infinite order. Some selected higher 
n-particle-71-hole correlations are also included. The in- 
clusion of triples correlations, that is the admixture of 
three-particle-three-hole correlations on the ground state 
of the nucleus with A nucleons and 3h — 2p excitations 
in the PR-EOM-CCM calculations of the A - 1 nucleus, 
on spectroscopic factors remains to be investigated. The 
largest uncertainty in our calculations is however most 
likely the effect of three-body interactions arising in chi- 
ral perturbation theory. A recent analysis by Otsuka 
et al [Slj demonstrates in shell-model calculations con- 
strained by the degrees of freedom of the sd-shell, that 
three-body interactions are important in order to obtain 
the experimental trend in binding energies for the oxy- 
gen isotopes. Three-body interactions, included as den- 
sity dependent corrections to the two-body interactions 
for the sd-shell, result in ^^O as unbound with respect to 
^^O. This is however still an unsettled topic. A similar 
effect can be obtained at the level of two-body effective 
interactions by including higher-lying single-particle ex- 
citations in many-body perturbation theory. The role of 
such three-body interactions in our calculations of spec- 
troscopic factors and single-particle energies is a topic for 
future investigations. 

Conclusion and Outlook We have computed single- 
particle energies and spectroscopic factors for hole states 



in ^ O using coupled-cluster theory at the level of singles 
and doubles correlations. The role of continuum states 
has also been included in our investigations. For the hole 
states the major influence of the continuum states is to 
give final single-particle energies and spectroscopic fac- 
tors which are almost independent of the chosen oscil- 
lator energy. The spectroscopic factors for protons and 
neutrons obtained with ab initio coupled-cluster calcula- 
tions support the emerging consensus that ^^0 is a dou- 
bly magic nucleus. In future work we plan to investigate 
the application of ab initio radial overlap functions to 
neutron knock-out reactions on drip line nuclei, thereby 
removing a level of model dependence from reaction stud- 
ies. 
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